G3: Genes, Genomes, Genetics

The Candida Genome Database (CGD; www.candidagenome.org) is both a model organism database and a fungal pathogen database. As a model organism database, CGD stores data for Candida albicans, which serves as a model species both for other Candida spp. and for non-Candida fungi that form biofilms and undergo routine morphogenic switching. As a fungal pathogen database, CGD now hosts locus pages for six species of the best-studied pathogenic fungi in the Candida group. Pathogenic Candida species have become increasingly drug resistant and there is thus a pressing need for research into basic Candida biology, epidemiology, phylogeny, and potential new antifungals, as well as a single location where all of the available data are collected, curated, and made easily searchable. CGD curates the gene-based Candida experimental literature in real time, extracting, organizing and standardizing gene annotations. CGD also links clinical data on disease to relevant Literature Topics to improve searchability for clinical researchers. Because CGD curates the literature for multiple species and most research focuses on aspects related to pathogenicity, we focus our curation efforts on assigning Literature Topic tags, collecting detailed mutant phenotype data, and assigning controlled Gene Ontology terms with accompanying evidence codes. Our Summary pages for each locus include the primary name and all aliases for that locus, a description of the gene and/or gene product, detailed ortholog information with links, a synteny view, a JBrowse window with a visual view of the gene on its chromosome, links to Phenotype, Gene Ontology, Interactions, and Expression pages, as well as sequence information, references cited on the summary page itself, and any locus notes. The database also serves as a community hub, where we link to various types of reference material of relevance to Candida researchers, including colleague information, news, and notice of upcoming meetings. We routinely survey the community to learn how the field is evolving and how needs may have changed. Here we describe CGDs new modern web interface and multiple new tools that have been added in the last 6 months, allowing, among other things, users to better understand the available expression data for a locus and seamlessly switch between species for a given locus.

The rat tapeworm, Hymenolepis diminuta, is an important laboratory model for uncovering molecular processes that underly the success of tapeworms as parasites. Despite its importance, a high-quality reference genome for this species is lacking. Here we present a highly contiguous and effectively complete genome of H. diminuta assembled from PacBio HiFi long-read sequencing data. Our primary assembly consists of 7 scaffolds (N50=29.25 Mb) with total length of 186.53 Mb, has only 7 gaps, and contains 95.7% complete Lophotrochozoan BUSCOs. Our assembly allows us to confirm aspects of Hymenolepis genome organization, such as high repeat content and unusual chromosomal ends, and to show that Hymenolepis genomes encode [~]10,000 genes. Together with annotations of nuclear tRNAs, mtDNA protein coding genes, and mtDNA tRNAs, our assembly currently provides one of the most complete genome resources for a tapeworm species and will enable research on parasitism, animal regeneration, development, and evolution.

The Tyr::CreERT2 transgenic mouse lines are essential tools for conditional gene manipulation in melanocytes and are widely used in melanoma research. Two independent lines are in common use: the Bosenberg line and the Larue line. Precise knowledge of transgene integration sites is critical for designing complex genetic crosses, yet the integration sites for these lines have only recently been characterised by whole genome sequencing. Here we report that a Tyr::CreERT2 mouse stock routinely used in BrafCA;Ptenflox melanoma models carries a previously undescribed integration on Chromosome 1, distinct from the previously reported Chromosome 2 integration. Using long-read nanopore sequencing, we mapped the integration to an intergenic locus between Alppl2 and Alpi, revealing a 2,430 bp genomic deletion at the insertion site. We developed position-specific junction PCR and qPCR assays to genotype this allele and confirmed that offspring are born at Mendelian ratios. Given that this stock is associated with elevated spontaneous melanoma penetrance, accurate characterisation of this allele has direct implications for the many laboratories employing this widely distributed melanoma model. SIGNIFICANCEThis study identifies a previously undescribed Tyr::CreERT2 transgene integration on Chromosome 1 in JAX strain 013590, a line widely used across the melanoma research community. Accurate characterisation of this allele has direct practical implications for the many laboratories that rely on this model for studies of melanoma initiation and progression.

Advances and accessibility to next-generation technology provide opportunities to sequence non-model organisms. Despite this increase in whole-genome sequencing data, annotation and the production of reference genomes remain limited. Reference genomes are a critical tool for a variety of studies in evolutionary biology, functional genomics, and conservation genetics. Tiger beetles (Cicindelidae) are a diverse and globally distributed family of beetles that serve as bioindicator taxa and flagship species for insect conservation. Here, we report highly complete, contiguous, and annotated genome assemblies representing draft reference genomes for three species of tiger beetle spanning the phylogeny. These draft reference genomes are for Audouins night-stalking tiger beetle, Omus audouini; the montane giant tiger beetle, Amblycheila baroni; and the western red-bellied tiger beetle, Cicindelidia sedecimpunctata. Article SummaryTiger beetles are a charismatic group with [~]3000 species distributed globally. Despite their popularity among insect enthusiasts and their role as bioindicators of ecosystem health, the group currently lacks a reference genome. This article outlines genome assembly and annotation for three tiger beetle species that span evolutionary relationships within the lineage. Quality control analyses show that the assemblies are reference quality and demonstrate high contiguity, completeness, and accuracy. The resulting draft genome annotations will be a valuable resource for scientific endeavors and allow for continued research on tiger beetles and their allies.

Gene model for the ortholog of tango (tgo) in the Sep. 2015 (UC Berkeley ASM127793v1/DbusGB1) Genome Assembly (GenBank Accession: GCA_001277935.1) of Drosophila busckii. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

We present a chromosome-level genome assembly of a female great grey owl (Strix nebulosa lapponica). The assembly comprises two pseudo-haplotypes of 1554 Mb and 1242 Mb, with 83.2% and 91.4% scaffolded into 40 autosomal chromosomes, in addition to the W and Z sex chromosomes both placed in hap1. Assembly completeness is high (BUSCO 99.2% and 94.8%), with 18,493 and 17,279 annotated protein-coding genes for hap1 and hap2, respectively. This genome establishes a reference for investigating genetic variation and chromosome evolution in great grey owls. Compared with the previous S. nebulosa assembly, this assembly includes both sex chromosomes, separates regions that were previously collapsed, and resolves 82 chromosomes total. While larger chromosomes show broadly conserved synteny across owl assemblies, the recovery of additional conserved microchromosome-associated genes suggests that ONT reads improved resolution of the smallest chromosomes relative to HiFi-based assemblies.

The potato leafhopper, Empoasca fabae (Harris, 1841), is a highly polyphagous, migratory insect pest of eastern North America that feeds on more than 200 herbaceous and woody plant species, causing substantial losses to forage and field crops. Despite its agricultural and ecological importance, no genome has been available for this species. Here, we present the first chromosome-level genome assembly of E. fabae, generated from Oxford Nanopore long reads, Illumina short reads, and Omni-C proximity-ligation data. The final assembly spans 908 Mb across 132 scaffolds, with 99.8% of the assembly captured in ten chromosome-length scaffolds (nine autosomes and an X chromosome) with a scaffold N50 of 96.2 Mb. The assembly is highly complete, recovering 92.4% of conserved hemipteran single-copy orthologs, and is composed of 47.6% repetitive sequence, dominated by long terminal repeat retrotransposons and unclassified elements. Read-depth comparison between male and female individuals supports assignment of a single sex-linked chromosome, consistent with an XO sex-determination system. BRAKER3 gene annotation predicted 31,406 protein-coding genes after retaining the longest isoform per locus. Comparative genome analysis against the two closest related Typhlocybinae species with genomes available, Matsumurasca onukii and Hebata decipiens, revealed extensive chromosome-scale collinearity, while defining a shared core gene repertoire. This reference genome provides a foundation for comparative and population genomic studies and for investigating genetic traits in this economically important crop pest species.

Multiple displacement amplification (MDA) enables whole-genome amplification from single cells, but introduces chimeric artifacts that severely compromise downstream analyses, particularly with long-read sequencing. Here, we systematically evaluate long-read PacBio HiFi sequencing of MDA amplified DNA from single cells using the model green alga Chlamydomonas reinhardtii. We show that MDA-derived libraries exhibit highly uneven coverage and extreme chimera rates impacting up to 70% of reads, leading to thousands of artefactual structural variants and misassemblies when assembled using algorithms designed for bulk sequencing. To overcome these challenges, we developed lrSAGA (long-read Single Amplified Genome Assembly), a novel tool to assemble long-read MDA sequencing datasets. Assemblies generated using lrSAGA are more complete, more contiguous, and have 75-95% fewer misassemblies compared to conventional assembly algorithms. Although overall contiguity is limited by MDA coverage dropouts, we demonstrate that up to 68% of the C. reinhardtii genome can be accurately assembled from just a single haploid cell. We further validated lrSAGA using published Oxford Nanopore and PacBio HiFi data from single or half Caenorhabditis elegans worms, generating accurate and highly complete assemblies. Applying our approach to single protist cells isolated from environmental water samples, we performed PacBio HiFi single-cell genome sequencing of four uncultivated microbial eukaryotes: an amoeboflagellate from the Naegleria genus, a flagellate from the Bodo genus, and two deep-branching flagellates from the enigmatic CRuMs supergroup, Collodictyon triciliatum and Diphylleia rotans. From single cells, we generated high-quality draft genome assemblies estimated to be 70-84% complete, demonstrating the potential of long-read single-cell genomics to unlock genome diversity from uncultivated microbial eukaryotes.

A highly contiguous and complete reference genome of Cercospora cf. flagellaris, the causal agent of foliar disease on many plant hosts including Cercospora leaf blight of soybean, was assembled using a combination of PacBio and Illumina sequencing reads. The genome assembly is 33.72 Mb in length and consists of 14 nuclear scaffolds and one mitochondrial contig. Four scaffolds have telomeric repeats on both ends and represent fully assembled chromosomes, while nine scaffolds represent partially assembled chromosomes with telomeric repeats on one end. The assembly has an N50 of 2.90 Mb and an L50 of 5 scaffolds. Genome annotation identified 11,268 genes, of which 947 and 360 were predicted to encode secreted proteins and effectors, respectively. Additionally, 512 genes were predicted to encode carbohydrate-active enzymes and 60 biosynthetic gene clusters were annotated. Taken together, this annotated genome assembly will be a valuable resource for genomics, host-pathogen interactions, and population biology research in this economically important pathosystem.

Gene model for the ortholog of Density regulated protein (DENR) in the Apr. 2013 (BCM-HGSC/Deug_2.0) (DeugGB2) Genome Assembly (GenBank Accession: GCA_000236325.2) of D. eugracilis. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Import of double-stranded RNA (dsRNA) into the germ line can have consequences that last for many generations. However, the role of such transgenerational regulation by extracellular dsRNA is unclear. In the nematode C. elegans, entry of dsRNA into the cytosol requires the transmembrane protein SID-1 and loss of SID-1 for a few generations causes changes in gene expression that can persist for hundreds of generations. Here we report an expanded number of such SID-1-dependent genes (SDGs) and analyze two germline-expressed SDGs: sdg-1 and sdg-2. Deleting sdg-1 reduces brood size in some lineages. An endogenous SDG-1::mCherry fusion protein shows conditional enrichment within nuclei, colocalization with perinuclear germ granules, and colocalization with microtubules. Although animals with SDG-1::mCherry have a normal brood size, they have fewer early progeny with some animals showing defective germline morphology. Deleting the sdg-1 open reading frame eliminates defects in most but not all the animals that express mCherry in a now sdg-1(-) background, suggesting transgenerational consequences of SDG-1::mCherry that persist in some siblings lacking sdg-1. Deleting sdg-2 also reduces brood size in some lineages. An endogenous SDG-2::mCherry fusion protein is constitutively detectable in the cytoplasm and nucleus. The sequence and predicted structure of SDG-2 suggest that it can interact with the Gli-type transcription factor TRA-1, which regulates spermatogenesis. Together, these results suggest that changes in SDG-1 or SDG-2 can impact reproduction. Therefore, the import of extracellular dsRNA or other SID-1 function(s) that regulate SDGs could have evolved to modulate the lingering impacts of ancestral epigenetic changes.

Gene model for the ortholog of raptor in the D. grimshawi May 2011 (Agencourt dgri_caf1/DgriCAF1) Genome Assembly (GenBank Accession: GCA_000005155.1) of Drosophila grimshawi. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of raptor in the D. eugracilis Apr. 2013 (BCM-HGSC/Deug_2.0) (DeugGB2) Genome Assembly (GenBank Accession: GCA_000236325.2) of Drosophila eugracilis. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Density regulated protein (DENR) in the May 2011 (Agencourt dgri_caf1/DgriCAF1) Genome Assembly (GenBank Accession: GCA_000005155.1) of D. grimshawi. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Genetic-background studies require defined perturbations that can be crossed reproducibly into many recipient backgrounds. We generated a Drosophila dilp2GS-rpr donor line for adult-inducible ablation of insulin-producing cells (IPCs), which secrete insulin-like peptides and provide a tractable model of insulin-deficient metabolic physiology. This line carries dilp2-GeneSwitch-GAL4 and UAS-reaper in cis on the same second chromosome homolog over a balancer. PCR genotyping and sequencing confirmed both transgenic elements in the candidate recombinant line. RU486 induction reduced dilp2 mRNA expression, supporting partial IPC ablation. Treatment-duration testing identified 8 days of RU486 as sufficient to increase whole-body glucose in the dilp2GS-rpr line but not in the background-matched control; food intake did not differ between RU486- and vehicle-treated flies. Across metabolic assays, whole-body glucose showed the clearest RU486- and line-dependent phenotype. This validated dilp2GS-rpr line enables testing how recipient genetic backgrounds modify inducible IPC/DILP metabolic phenotypes and provides a framework for similar linked donor-line resources.

The coffee bean weevil, Araecerus fasciculatus (Coleoptera, Curculionoidea, Anthribidae), is a cosmopolitan pest of over 100 stored agricultural commodities, with particular economic impact on coffee (Coffea arabica). Although two chromosome-level anthribid genomes have recently been released as part of the Darwin Tree of Life (DToL) project (Booth et al. 2024; Crowley et al. 2025), no functionally annotated genome has been available for the family. Here we present a draft genome assembly for A. fasciculatus, generated from PacBio HiFi long reads and processed through a three tiered metagenomic filtering pipeline to remove host plant (C. arabica) and microbial contamination. The final assembly spans 475 Mb across 3,617 scaffolds (N50 = 170 kb) with 88.5% BUSCO completeness (insecta_odb10) and only 3.1% duplication. Gene prediction with BRAKER2 identified 22,384 protein-coding genes, of which 11,783 received functional annotations through SwissProt similarity. Notably, we identified 92 cytochrome P450 (CYP) genes, including tandem gene clusters on two scaffolds (4 genes on ptg000464l, 5 genes on ptg001867l), suggestive of lineage-specific expansion through tandem duplication. Homology searches against Drosophila melanogaster caffeine-metabolizing P450s (CYP12D1, CYP6d5, CYP6a8) recovered strong matches (e-values 9.7 x 10-110 to 5.4 x 10-101, 33-38% identity). In stark contrast, comprehensive BLAST searches for bacterial caffeine N-demethylase genes (ndmA/B/C/D), which mediate caffeine degradation via horizontal gene transfer in the coffee berry borer Hypothenemus hampei (Scolytinae), returned zero hits across the A. fasciculatus genome, predicted proteome, and associated bacterial scaffolds. AlphaFold2 structure prediction of four top Araecerus P450 candidates produced high-confidence models (pLDDT 84.5-93.9, pTM 0.735-0.930) with conserved P450 catalytic motifs. Foldseek structural homology searches confirmed that all four candidates adopt cytochrome P450 folds (top hits: human CYP3A4, CYP3A7, CYP11A1; TM-scores 0.90-0.92; probability 1.000), with zero hits to bacterial Rieske-fold enzymes. Molecular docking of caffeine against these structures yielded binding affinities of -5.41 to -5.80 kcal/mol for the Araecerus candidates, comparable to or exceeding the -5.55 kcal/mol obtained for the experimentally validated Drosophila CYP6a8 and substantially stronger than the -3.70 kcal/mol for the bacterial NdmA structural outgroup (PDB: 6ICP). Phylogenetic analysis revealed that all four candidates have clear orthologs in two non-seed-feeding DToL anthribids (Pseudeuparius sepicola and Platystomos albinus), demonstrating that these P450 genes predate the dietary transition to caffeine-containing seeds. The Araecerus candidates predominantly belong to the CYP6 family (clan 3), whereas the primary Drosophila caffeine P450 CYP12D1 belongs to the mitochondrial clan, confirming convergent recruitment of different P450 subfamilies for caffeine metabolism. These results support the hypothesis that A. fasciculatus employs an insect-encoded, P450-mediated caffeine detoxification pathway fundamentally distinct from the bacterial horizontal gene transfer mechanism documented in Scolytinae. This represents convergent evolution of caffeine resistance via independent molecular strategies within Curculionoidea, and provides the first functionally annotated genomic resource for comparative studies across the Anthribidae.

Deletion of the Drosophila-specific gene Neprilysin-like 15 (Nepl15) results in significant reductions in glycogen and glycerolipid storage in mutant males and increased glycogen storage in mutant females without affecting food intake. Previous studies also indicated downregulation of insulin/mTOR signaling, the central pathway regulating nutrient homeostasis, growth, fertility, locomotor activity, and lifespan. Consistent with these metabolic alterations, Nepl15 mutants exhibit several anti-obesity and healthy-aging phenotypes but display markedly reduced survival under starvation, suggesting impaired nutrient reserve utilization. The present study investigated the intracellular mechanisms underlying these phenotypes by examining insulin signaling and the expression of genes involved in carbohydrate and lipid metabolism. Although transcript levels of the insulin-like peptides (Dilp2, Dilp5, and Dilp6) and the insulin receptor (InR) remained unchanged, Nepl15 mutants exhibited reduced Akt phosphorylation, increased dFoxo abundance, reduced Glut1 expression, and previously reported suppression of mTOR signaling, indicating attenuation of insulin signaling downstream of the insulin receptor. Consistent with impaired anabolic signaling, male mutants showed reduced expression of glycogen metabolic genes (GlyS and GlyP), whereas genes involved in lipid metabolism exhibited predominantly male-specific reductions, including Lipin, Acc, Fasn1, and Fasn2, while Midway and Brummer remained unchanged. In addition, expression of the metabolic regulator PGC1 (spargel) was significantly reduced in males, consistent with previously reported reductions in AMPK expression. Together, these findings demonstrate that Nepl15 functions as an upstream regulator of insulin-dependent anabolic signaling and nutrient partitioning, linking intracellular signaling to glycogen and lipid storage independently of nutrient intake. These results identify Nepl15 as a previously unrecognized regulator of metabolic homeostasis in Drosophila and provide new insight into the metabolic functions of the evolutionarily conserved neprilysin family.

Gene model for the ortholog of tango (tgo) in the May 2011 (Broad dper_caf1/DperCAF1) Genome Assembly (GenBank Accession: GCA_000005195.1) of Drosophila persimilis. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

The cultivated dahlia has a complex evolutionary history that is likely the driver for its diverse floral forms. Previous studies have proposed cultivated dahlia as an allotetraploid, autotetraploid, and in one study an auto-allooctoploid. Here, we used PacBio HiFi long-reads and Dovetail Omni-C reads to assemble and annotate a genome for a tetraploid Dahlia variabilis cultivar Edna C. We identified Edna C as a segmental tetraploid with signature of auto- and allo-polyploidy events based on subgenome-specific repeat partitioning and single-copy gene tree construction within homoeologous groups. Additionally, we identified a cultivated dahlia-specific centromeric monomer shared between Dahlia variabilis Edna C and the Dahlia pinnata Kelvin Floodlight. Based on finding three major chromosome clusters, we propose that Dahlia variabilis Edna C is a segmental tetraploid with evidence of auto- and allopolyploidy.

During meiosis, accurate chromosome segregation requires significant condensation and compaction. These processes are mediated by condensins, cohesins, and histone tail modifications. We identified that MET-2, a histone methyltransferase that catalyzes the dimethylation of histone H3 lysine 9 (H3K9me2), differentially impacts chromosome size in the male vs. female C. elegans germline. In met-2 null worms, autosomes during spermatogenesis are significantly larger than wild type, while chromosome size during oogenesis is unaffected. X-univalent size in males is also unaffected by loss of MET-2, indicating MET-2 differentially regulates autosomal and X-chromosome compaction in male spermatogenesis. Autosome size is not changed when males harbor a catalytically deficient MET-2 (met-2CD) or have mutations preventing germline histone H3K9 methylation (H3K9R). In addition, met-2 males, in contrast to met-2CD or H3K9R males, have more active RNA pol II in later stages of meiosis. These data suggest MET-2 plays a noncatalytic role in mediating chromosome structure and transcription. In met-2 male germ lines, genes on the X chromosome, which is typically enriched in H3K9me2, are significantly more likely to be upregulated than genes on autosomes, even though X-univalent size is unchanged. These results suggest that MET-2 plays a sex-specific role that is not limited to its enzymatic activity.